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15 Abstract 

16 A GEOS V 2 CCM simulation with an internally generated quasi-biennial oscillation (QBO) 

17 signal is compared to an otherwise identical simulation without a QBO. In a present-day 

18 climate, inclusion of the modeled QBO makes a significant difference to stratospheric dynamics 

19 and ozone throughout the year. The QBO enhances variability in the tropics, as expected, but 

20 also in the polar stratosphere in some seasons. The modeled QBO also affects the mean 

21 stratospheric climate. Because tropical zonal winds in the baseline simulation are generally 
easterly, there is a relative increase in zonal wind magnitudes in tropical lower and middle 


22 
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23 stratosphere in the QBO simulation. Extra-tropical differences between the QBO and ‘no QBO’ 

24 simulations thus reflect a bias toward the westerly phase of the QBO: a relative strengthening 

25 and poleward shifting the polar stratospheric jets, and a reduction in Arctic lower stratospheric 

26 ozone. 

27 
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28 1 Introduction 

29 The quasi-biennial oscillation (QBO) is the leading mode of variability in the tropical lower and 

30 middle stratosphere [Baldwin et al., 2001]. The QBO is characterized by a downward- 

31 propagating pattern of alternating easterly and westerly zonal winds in the equatorial region, with 

32 a period of 28 months (Figure la), and is driven by gravity and planetary-scale waves. The 

33 zonal wind QBO induces changes in the tropical stratospheric circulation, affecting the 

34 concentrations of ozone and other trace constituents [Gray et al., 1989; Butchart et al., 2003; 

35 Tian et al., 2006]. 

36 

37 The phase of the QBO affects the polar stratosphere. Holton and Tan [1980] and Lu et al. [2008] 

38 showed that the phase of the QBO modulates the strength of the Arctic vortex in mid-winter; the 

39 vortex is weakest during the easterly phase of the QBO. Hurwitz et al. [2011a, 201 lb] found that 

40 the phase of the QBO modulates the strength of the Antarctic vortex in austral summer, during 

41 warm pool El Nino (WPEN) events [see Kug et al., 2009], Analogously to the Arctic response, 

42 the Antarctic vortex is weakest during WPEN events coincident with the easterly phase of the 

43 QBO. In addition, QBO phase-related differences in the strength of the polar vortices modulate 

44 polar ozone loss [Lait, 1989; Randel and Cobb, 1994]. However, as the QBO signal is intrinsic 

45 to the observational record, the time-averaged impact of the easterly and westerly phases of the 

46 QBO on mean zonal wind, temperature and ozone cannot be evaluated using atmospheric data. 

47 

48 A chemistry-climate model (CCM) is the ideal tool for understanding the net impact of the QBO. 

49 Punge and Giorgetta [2008] quantified the net effect of the QBO on the late 20 th century 

50 stratospheric climate by comparing two CCM simulations: one without a QBO signal, and the 
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51 other with zonal winds between 90 hPa and 10 hPa nudged to profiles taken at Singapore. The 

52 authors found that inclusion of the QBO signal made a significant difference to zonal winds, 

53 temperature and ozone, mainly in the deep tropics. The vertical pattern of changes in tropical 

54 upwelling was consistent with changes in the transport of trace species. However, the value of 

55 the conclusions reached by Punge and Giorgetta [2008] is limited because the authors tested a 

56 prescribed QBO. Nudging to observed zonal winds does not allow the tropical and mid-latitude 

57 stratosphere to adjust realistically to changes in e.g., the QBO phase, nor does it allow for full 

58 interaction between stratospheric ozone and climate. 

59 

60 The Goddard Earth Observing System Chemistry-Climate Model, Version 2 (GEOS V2 CCM) 

61 can be run with or without an internally-generated QBO. In the formulation of the GEOS V2 

62 CCM evaluated by SPARC CCMVal [2010], both tropical lower stratospheric variability and the 

63 QBO amplitude were negligible, typical of CCMs that lack a nudged or internal QBO signal. In 

64 contrast, a more recent model formulation [introduced by Hurwitz et al., 2011b] has an 

65 internally-generated QBO with a realistic periodicity and amplitude. Comparing the two 

66 formulations of the GEOS V2 CCM, this paper quantifies the net effect of an internally 

67 generated QBO on stratospheric climate and variability. Section 2 provides a brief description of 

68 the two formulations of the CCM, as well as the simulations used to test the net impact of an 

69 internal QBO on the stratosphere. The effects of the internal QBO on the mean and variance of 

70 zonal winds, temperature and ozone are shown in Section 3. Section 4 provides a summary and 

7 1 discussion of the implication of the results for future CCM studies. 

72 

2 


73 


Model and simulations 
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74 This paper considers the net impact of the QBO in the GEOS V2 COM. The GEOS V2 CCM 

75 couples the GEOS-5 general circulation model (GCM) [Molod et al., 2011] with a 

76 comprehensive stratospheric chemistry module [Pawson et al., 2008]. The model has 2° latitude 

77 x 2.5° longitude horizontal resolution and 72 vertical layers, with a model top at 0.01 hPa. 

78 Predicted distributions of water vapor, ozone, greenhouse gases (CO2, CH4, and N2O) and CFCs 

79 (CFC-11 and CFC-12) feedback to the radiative calculations. The present formulation of the 

80 GEOS V2 CCM is the same as in Hurwitz et al. [2011b]. The GEOS V2 CCM performed well 

81 in the SPARC CCMVal [2010] detailed evaluation of stratospheric processes. 

82 

83 The model’s gravity wave parameterization computes the momentum and heat deposition to 

84 orographic and non-orographic gravity wave-breaking using the linear saturation theory by 

85 Lindzen [1981], Orographic gravity wave stress is specified using the formulation derived by 

86 McFarlanc [1987] and given at the top of the subgrid-scale mountains. Subgrid-scale orography 

87 is assumed to be horizontally isotropic. Heat transfer due to gravity wave-breaking is computed 

88 from the deposition of the gravity wave energy flux into the mean flow following Warner and 

89 McIntyre [2001] and Shaw and Shepherd [2009]. In addition, for the conservation of angular 

90 momentum and energy, gravity waves stress and energy flux is gradually dissipated in the top 

91 five model layers, as suggested by Shaw and Shepherd [2007]. 

92 

93 Modeled variability in tropical stratospheric zonal winds depends on the details of the non- 

94 orographic gravity wave drag (GWD) scheme. Early versions of the GEOS V2 CCM did not 

95 generate a QBO i.e., zonal winds in the equatorial lower stratosphere were generally easterly [see 
Figures lb and 2e-h; SPARC CCMVal, 2010]. A new formulation of the non-orographic 


96 



97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 


M. M. Hurwitz et al. (2011) 
Stratospheric Influence of a Modeled QBO 

6 

gravity wave drag scheme allows the model to generate a spontaneous QBO with a realistic 
period (30 months at 30 hPa and 50 hPa; see Figure lc) and zonal wind amplitude [see Hurwitz 
et al., 201 lb]. As non-orographic gravity waves often accompany precipitation (e.g., convective 
and frontal systems; see Richter et al. [2010]), the latitudinal structure of the gravity wave 
spectrum is designed to mimic the structure of the climatological mean precipitation field. A 
tropical peak in non-orographic gravity wave stress is necessary for the generation of an internal 
QBO in the model. A 700-km wavelength is used for the tropical non-orographic waves to 
prevent an excessive downward propagation of the semi-annual oscillation into the lower 
stratosphere, and thus contamination of the QBO signal. 

Two time-slice simulations will be considered in Section 3 of this paper: One 50-year 

simulation with an internal QBO (hereafter “QBO”) and one 20-year simulation without a QBO 
(“NQ”). Both simulations are forced by set of SST and sea ice climatologies, each with a 
repeating annual cycle, with conditions composited from 10 ENSO neutral (NTRL) years that 
span the satellite era [as in Hurwitz et al., 201 lb]. HadlSSTl SSTs and sea ice concentrations at 
1° x 1° resolution [Rayner et al., 2003] are used to prepare the composites. All simulations used 
fixed greenhouse gas and ozone-depleting substance boundary conditions representative of the 
year 2005. Variability related to the solar cycle and volcanic eruptions is not considered. 

Modeled temperature, zonal winds and ozone are compared with a composite of 12 observed 
NTRL years distributed throughout the 1979-2010 period, using the Modem Era Retrospective- 
Analysis for Research and Applications (MERRA). MERRA is a reanalysis dataset based on an 
extensive set of satellite observations and on the Goddard Earth Observing System Data Analysis 
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System, Version 5 (GEOS-5) [Bosilovich et al., 2008; Rienecker et al., 2011], The MERRA 
reanalysis has vertical coverage up to 0.1 hPa, and for this study, is interpolated to 1.25° x 1.25° 
horizontal resolution. 

3 Results 

3.1 QBO influence on zonal winds and temperature 

Seasonal mean zonal wind climatologies in the NQ simulation are shown in Figures 2e-h. 
Above 100 hPa, tropical zonal winds are easterly throughout the year in this simulation. The 
tropical zonal wind magnitude maximizes in the upper stratosphere in the DJF and JJA seasons. 
The stratospheric polar night jets maximize in the winter season, in the respective hemisphere. 
The modeled winds compare well with the MERRA reanalysis (Figures 2a— d); in NQ, the polar 
jets are weaker than observed while tropical stratospheric easterlies are stronger than observed 
(Figures 2e-h). 

Inclusion of the modeled QBO makes a significant difference to zonal winds in the stratosphere. 
Figures 2i-l show zonal wind differences between the QBO and NQ simulations. By definition, 
in the deep tropics, the modeled QBO increases mean zonal wind magnitudes throughout the 
year, with the largest increases of 10-20 ms" 1 in the middle and upper stratosphere. At high 
latitudes, the sign of the zonal wind differences is generally independent of height: negative 
differences equatorward of 60° and positive differences poleward of 60° indicate poleward 
shifting of the polar stratospheric jets. Poleward shifting of the Southern Hemisphere jet is 
evident throughout the year, while in the Northern Hemisphere it is limited to the transition 
seasons. During the DJF season, the Arctic jet strengthens by 5 m s’ 1 (in the lower stratosphere) 
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143 to 15 m s' 1 (in the upper stratosphere), consistent with the increased frequency of QBO-westerly 

144 years in the QBO simulation. 

145 

146 The modeled QBO affects stratospheric zonal wind variability. Zonal wind variability increases 

147 significantly in the QBO simulation, as compared with NQ, throughout the tropical stratosphere 

148 and in all seasons (Figures 2m-p). The variance of tropical zonal winds increases tenfold in the 

149 lower stratosphere, and by 3-5 times in the middle and upper stratosphere. Decreased mid- 

150 latitude variability is evident in the DJF and JJA seasons in the SH and in the SON season in the 

151 NH. 

152 

153 The GEOS V2 CCM reproduces the Holton-Tan (1980) relation between the phase of the QBO 

154 and Arctic vortex strength in winter. In the QBO simulation, 30 hPa zonal winds in the 

155 equatorial region are strongly positively correlated with zonal winds in the Arctic stratosphere in 

156 November, December and January: the modeled Arctic vortex is relatively weaker during the 

157 easterly phase of the QBO (not shown). Furthermore, the modeled QBO enhances polar 

158 variability in zonal winds, particularly during the spring and summer seasons. Enhanced 

159 variability in the strength of the polar vortices in late winter and spring implies more variability 

160 in the timing of the vortex breakup, and thus, more variability in polar ozone. Links between 

161 vortex strength and ozone will be discussed further in the next section. 

162 

163 ENSO neutral stratospheric temperatures are well simulated by the model (Figures 3a-h). 

164 Consistent with the weaker-than-observed zonal winds (Figures 2a-h), temperatures in the 

165 extra-tropical winter stratosphere are biased high in the NQ simulation (i.e., the polar vortices 
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166 are weaker than in MERRA). The modeled QBO acts to strengthen the polar vortices, cooling 

167 the extra-tropical stratosphere in autumn and winter up to 10 K (Figures 3i-l). In addition, 

168 decreased temperatures are seen in the Arctic lower stratosphere in boreal spring (Figure 3j). 

169 Inclusion of the QBO increases temperatures on the equatorward side of the polar jets, consistent 

170 with the poleward shifting of the polar vortices. 

171 

172 QBO-related changes in temperature variability maximize in the middle stratosphere (Figures 

173 3m-p). Tropical variability increases in all seasons. Variability in the polar stratosphere 

174 generally increases as a result of the QBO, while there is a region of reduced variability around 

175 60°S, associated with the stronger and more stable Antarctic vortex. 

176 

177 3.2 QBO influence on ozone 

178 Seasonal mean ozone mixing ratio fields, in the NQ simulation, are shown in Figures 4e-h. 

179 Ozone maximizes at approximately 10 ppmv in the deep tropics, at 10 hPa, and decreases with 

180 latitude. There is little inter-seasonal variability in tropical ozone. Modeled ozone is generally 

181 consistent with the MERRA reanalysis, though modeled peak ozone values in the equatorial 

182 middle stratosphere are up to 1 ppmv lower than in MERRA (Figures 4a-h). The MERRA 

183 ozone product is generally in agreement with MLS v3.3 ozone [Froidevaux et al., 2008; Jiang et 

184 al., 2007; Livesey et al., 2008] during the 2005-2010 period, with a slight high bias in the 

185 tropical middle stratosphere and low bias in the polar upper stratosphere. 

186 

187 QBO - NQ ozone differences are statistically significant throughout the stratosphere (Figures 4i- 

188 1). In all seasons, the peak in ozone mixing ratio around 10 hPa is reduced by ~1 ppmv in the 
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189 QBO simulation as compared with NQ. These ozone differences reflect the relative warming of 

190 the tropical middle stratosphere in the QBO simulation (see Figures 3i-l and Figure 5). In the 

191 Arctic lower stratosphere, inclusion of the QBO reduces ozone in the SON, DJF and MAM 

192 seasons; this result is consistent with the relative strengthening of the Arctic vortex and increase 

193 in chemical ozone depletion (as discussed in Section 3.1). During the SON or “ozone hole” 

194 season, inclusion of the modeled QBO increases ozone mixing ratios (i.e., reduces ozone 

195 depletion) in the Antarctic, reflecting increased temperatures in the polar lower stratosphere. In 

196 other seasons, the impact of the modeled QBO on Antarctic ozone is not as straightforward: the 

197 sign of the mean change in ozone varies as a function of latitude and altitude. 

198 

199 Inclusion of the modeled QBO has a significant impact on ozone variability (Figure 4m-p). 

200 Ozone variability increases in the tropics and subtropics in all seasons, and in the Arctic lower 

201 stratosphere in DJF, MAM and JJA. Decreased variability is restricted to the uppermost polar 

202 stratosphere in spring and summer. 

203 

204 4 Discussion 

205 This paper showed that a model’s representation of the QBO makes a significant difference to 

206 stratospheric climate and variability. Comparing a GEOS V2 CCM simulation with a realistic, 

207 internally generated QBO to one without, the addition of a QBO enhances tropical variability in 

208 zonal wind, temperature and ozone throughout the year. Extra-tropical zonal wind variability is 

209 also enhanced, mainly in the spring and summer seasons. 


210 
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211 Because the QBO is by definition an oscillating phenomenon, one might assume that a modeled 

212 QBO does not exert any net impact on stratospheric climate. This work tested and rejected this 

213 assumption: 

214 (1) Adding an internal QBO signal affects the mean stratospheric climate. 

215 In the GEOS V2 CCM, inclusion of a QBO signal shifts both polar vortices further poleward. 

216 The polar vortices are strengthened, particularly in autumn and winter. The QBO reduces 

217 ozone concentrations in the tropical middle stratosphere and Arctic lower stratosphere. 

218 (2) The mean response to the QBO depends on the background zonal wind field, in the model 

219 version without a QBO signal. 

220 Because tropical zonal winds in the NQ simulation are generally easterly, there is a relative 

221 increase in zonal wind magnitudes in tropical lower and middle stratosphere in the QBO 

222 simulation. Extra-tropical differences between the QBO and NQ simulations reflect a bias 

223 toward the westerly phase of the QBO: a relative strengthening and poleward shifting the 

224 polar stratospheric jets, and a reduction in Arctic lower stratospheric ozone. The annual 

225 mean impact of the QBO on the polar stratosphere is larger in the GEOS V2 CCM (up to 12 

226 m s' 1 , in both hemispheres) than in the MAECHAM4-CHEM CCM (no significant zonal 

227 wind differences) [Punge and Giorgetta, 2008], likely reflecting larger QBO - ‘no QBO’ 

228 tropical zonal wind differences in the GEOS V2 CCM and/or increased statistical robustness 

229 due to the greater length of the GEOS V2 CCM simulations. 

230 

231 At 10 hPa, equatorial ozone is negatively correlated with temperature (Figure 4). This result 

232 suggests that negative QBO - NQ ozone differences (e.g., Figure 4i) in this region can be 

233 explained by the relatively higher temperatures in the QBO simulation. Furthermore, the linear 
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234 relation between ozone and temperature in the tropical middle stratosphere suggests that, in a 

235 climate change simulation, inclusion of the QBO would likely affect the modeled representation 

236 of observed ozone values but would not affect ozone trends. 

237 
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337 Figure Captions 

338 

339 Figure 1: Zonal wind frequency spectra for 4°S-4°N in (a) MERRA, 1979-2009, (b) the NQ 

340 simulation, and (c) the QBO simulation. 

341 

342 Figure 2: Seasonal and zonal mean zonal wind, as a function of latitude and altitude, in the 

343 NTRL composite in MERRA (a, b, c, d; m s' 1 ) and in the NQ simulation (e, f, g, h; m s' 1 ), (i, j, 

344 k, 1; m s' 1 ) Seasonal mean zonal wind differences between the QBO and NQ simulations, as a 

345 function of latitude and altitude, (m, n, o, p) Ratio of the variance between the QBO and NQ 

346 simulations, as a function of latitude, altitude and season. White contours indicate no difference 

347 in the mean (variance). Black Xs indicate regions where differences in the mean (variance) are 

348 significant at the 95% confidence level. 

349 

350 Figure 3: Seasonal and zonal mean temperature, as a function of latitude and altitude, in the 

351 NTRL composite in MERRA (a, b, c, d; K) and in the NQ simulation (e, f, g, h; K). (i, j, k, 1; K) 

352 Seasonal mean zonal wind differences between the QBO and NQ simulations, as a function of 

353 latitude and altitude, (m, n, o, p) Ratio of the variance between the QBO and NQ simulations, as 

354 a function of latitude, altitude and season. White contours indicate no difference in the mean 

355 (variance). Black Xs indicate regions where differences in the mean (variance) are significant at 

356 the 95% confidence level. 

357 

358 Figure 4: Seasonal and zonal mean ozone mixing ratio, as a function of latitude and altitude, in 

359 the NTRL composite in MERRA (a, b, c, d; ppmv) and in the NQ simulation (e, f, g, h; ppmv). 
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360 (i, j, k, 1; ppmv) Seasonal mean zonal wind differences between the QBO and NQ simulations, as 

361 a function of latitude and altitude, (m, n, o, p) Ratio of the variance between the QBO and NQ 

362 simulations, as a function of latitude, altitude and season. White contours indicate no difference 

363 in the mean (variance). Black Xs indicate regions where differences in the mean (variance) are 

364 significant at the 95% confidence level. 

365 

366 Figure 5: Temperature versus ozone at 10°S-10°N, 10 hPa in January/February in the QBO 

367 simulation (triangles) and NQ simulation (squares). The thick black line denotes the linear fit 

368 between temperature and ozone in the QBO simulation. 
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Figure 1: Zonal wind frequency spectra for 4°S— 4°N in (a) MERRA, 1979-2009, (b) the NQ 
simulation, and (c) the QBO simulation. 
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Figure 2: Seasonal and zonal mean zonal wind, as a function of latitude and altitude, in the 
NTRL composite in MERRA (a, b, c, d; m s' 1 ) and in the NQ simulation (e, f, g, h; m s' 1 ), (i, j, 
k, 1; m s' 1 ) Seasonal mean zonal wind differences between the QBO and NQ simulations, as a 
function of latitude and altitude, (m, n, o, p) Ratio of the variance between the QBO and NQ 
simulations, as a function of latitude, altitude and season. White contours indicate no difference 
in the mean (variance). Black Xs indicate regions where differences in the mean (variance) are 


significant at the 95% confidence level. 
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Figure 3: Seasonal and zonal mean temperature, as a function of latitude and altitude, in the 
NTRL composite in MERRA (a, b, c, d; K) and in the NQ simulation (e, f, g, h; K). (i, j, k, 1; K) 
Seasonal mean zonal wind differences between the QBO and NQ simulations, as a function of 
latitude and altitude, (m, n, o, p) Ratio of the variance between the QBO and NQ simulations, as 
a function of latitude, altitude and season. White contours indicate no difference in the mean 
(variance). Black Xs indicate regions where differences in the mean (variance) are significant at 


the 95% confidence level. 
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Figure 4: Seasonal and zonal mean ozone mixing ratio, as a function of latitude and altitude, in 
the NTRL composite in MERRA (a, b, c, d; ppmv) and in the NQ simulation (e, f, g, h; ppmv). 
(i, j, k, l; ppmv) Seasonal mean zonal wind differences between the QBO and NQ simulations, as 
a function of latitude and altitude, (m, n, o, p) Ratio of the variance between the QBO and NQ 
simulations, as a function of latitude, altitude and season. White contours indicate no difference 
in the mean (variance). Black Xs indicate regions where differences in the mean (variance) are 


significant at the 95% confidence level. 
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Figure 5: Temperature versus ozone at 10°S-10°N, 10 hPa in January/February in the QBO 
simulation (triangles) and NQ simulation (squares). The thick black line denotes the linear fit 
between temperature and ozone in the QBO simulation. 



